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ABSTRACT 


The  transfer  of  kinetic  energy  by  individual  wave 
numbers  is  computed  at  various  latitudes  using  data  ob¬ 
tained  from  a  harmonic  analysis  of  the  hemispheric  500 
millibar  pressure  height.  Linear  correlations  are  sought 
between  the  computed  values  of  kinetic  energy  transfer  and 
computed  24-hour  changes  of  kinetic  energy  of  the  pertur¬ 
bations  for  the  same  and  for  adjacent  wave  numbers. 

The  writer  wishes  to  express  sincerest  appreciation 
to  Professor  Frank  L.  Martin  for  his  suggestion  of  the 
topic  and  for  his  guidance  during  the  study. 
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1. 


Introduction. 


One  phase  in  the  study  of  the  general  circulation  of 
the  atmosphere  has  been  the  investigation  of  mechanisms  by 
which  the  kinetic  energy  balance  is  maintained  in  spite  of 
continual  frictional  dissipation  Widger  [14] ,  Kao  [6] , 
and  others  have  established  that  a  mean  northward  trans¬ 
port  of  kinetic  energy  in  mid-latitudes  is  effected  by  the 
large  scale  perturbations  in  the  atmosphere  Saltzman 
[9,  10]  has  investigated  the  spectral  distribution  of  per¬ 
turbation  kinetic  energy.  Saltzman  and  Fleisher  [11]  have 
studied  the  rate  of  transfer  of  kinetic  energy  between  the 
perturbations  and  the  mean  flow  on  a  daily  basis,  while 
Benton  and  Kahn  [1]  have  conducted  a  similar  investigation 
on  a  seasonal  basis.  Saltzman  and  Fleisher  [12]  have  also 
estimated  the  interchange  of  eddy  kinetic  energy  between 
groups  of  perturbations,  assuming  that  the  total  kinetic 
energy  interchange  is  zero. 

It  is  the  purpose  of  this  study  to  correlate  the  hori¬ 
zontal  geostrophic  transfer  of  kinetic  energy  by  wave  num¬ 
ber  with  24-hour  changes  of  kinetic  energy  of  the  same  and 
different  wave  numbers  existing  in  the  same  and  nearby 
latitude  belts.  The  data  used  are  from  a  Fourier  analysis 
of  500-mb  contours  at  successive  5-degree  latitude  belts. 
This  Fourier  analysis,  which  covers  the  winter  season  of 
1947-48,  was  made  available  by  the  U.  S„  Navy  Weather 
Research  Facility,  Norfolk.  Further  discussion  concerning 
the  data  appears  in  a  later  section. 
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2.  The  energy  transfer  function. 

Using  the  equations  of  quasi-horizontal  motion, 
Saltzman  [9]  has  developed  an  expression  for  the  change  of 
eddy  kinetic  energy  within  a  closed  mass  of  air, M  ,  e.g. 
the  entire  atmosphere,  which  is 


r 


3 

at 


KE  s 


J 

tA 


V'  £2ii  ll-^\  -H  vV 


u  V 


a  So 


-  uV^tan4)  -4-u'gj'^  + 


•f  V  ^ 

c>-p 


1 — 7- 


Su) 


_  jc^\/'oYl'  ~  If*  s 

M  /A 


(2.1) 


Here  (f)  is  latitude,  7t  longitude,  ^  pressure,  and  Q  ,  the 
earth-radius.  The  eastward  and  northward  components  of  the 
wind  are  U  and  V  respectively,  while  the  vertical  motion 
is  oJ  .  Primes  are  associated  with  perturbation  quantities 
and  bars  with  averages  around  a  latitude  circle.  Thus,  for 
example , 

(  ^  -  in  I ^ 

o 

indicates  the  barred-type  average.  The  mass  element  in 
equation  (2.1)  will  be  discussed  in  a  later  section. 
Saltzman  also  shows  that  the  change  in  the  mean  kinetic 
energy,  ^  f  KE  '  contains  the  first  integral  on  the 

right  in  (2.1),  but  with  the  opposite  sign  affixed. 

It  may  be  seen  that  the  first  integral  in  equation 
(2.1)  contains  terms  which  are  products  of  the  macroscopic 
Reynolds  stresses  and  the  derivatives  of  the  mean  flow. 
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The  remaining  integrals  represent  contributions  to  the  eddy 
kinetic  energy  by  conversion  of  "eddy  available  potential 
energy"  and  frictional  dissipation,  respectively.  Inasmuch 
as  this  study  is  concerned  only  with  the  effects  of  hori¬ 
zontal  transfer  of  kinetic  energy,  no  further  discussion  of 
the  geopotential  or  frictional  terms  will  be  entertained. 

If  the  geostrophic  approximation  is  made,  as  is  done 
in  this  study,  the  mean  meridional  component  ~\/'  will  vanish. 
Furthermore,  since  it  is  impossible  to  determine  vertical 

motions  from  the  Fourier  analysis  of  the  height  field  at  a 

.  / 

single  level,  those  terms  involving  CJ  will  be  neglected. 

With  these  simplifications,  the  eddy  kinetic  energy 
change  may  be  expressed  in  the  form 


For  an  interpretation  of  the  integral  appearing  in 
equation  (2.2),  a  vertical  element  of  area,  Sz  -  acos<t>6'^h2 
contained  in  a  latitude  wall,  will  first  be  considered. 

Here  the  eastward  momentum  per  unit  mass  of  air  is 

u  =  u  +  u' 

Any  transport  of  eastward  momentum  normal  to  Ss  must  be 
brought  about  by  the  meridional  wind  component 

V  -  ^ 

With  the  geostrophic  approximation,  the  longitudinally- 
averaged  meridional  transport  of  eastward  momentum  due  to 
large  scale  turbulence  is  given,  after  Haltiner  and 
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Martin  [4] , 


LiV  =  a^){v^  -  A/^)  ^ 

The  integral  in  equation  (2.2)  will  be  denoted  T"  , 
and  will  be  termed  the  "eddy  transfer  function. "  It  will 
eventually  involve  contributions  to  the  kinetic  energy 
change  by  different  wave  numbers.  A  positive  value  for 
'~V  represents  an  exchange  of  kinetic  energy  from  the  mean 
flow  into  that  of  the  eddy  flow. 

3.  The  energy-transfer  spectrum. 

Since  U.  and  V  may  be  considered  real  and  continuous 
functions  of  longitude  A-  ,  each  having  unique  values  at  all 
points  of  a  given  latitude  circle,  they  may  be  expanded  in¬ 
to  Fourier  series  of  the  form 


U(A)  = 


ajo) 


cos  nTv  + 


sin  OA 


(3.1) 


where 


Likewise 


2.Tr 


211 


c 


/ 


^TT 


b„(h^  -  IT  I 


(3.2) 


V(7v)  = 


a  (0) 

AT 


+  7  a  (n)  cos  OA  + 


b^(n)  sin 


n7\ 


,0  0/ 


where  ,  Clv^'n)  and  are  defined  similarly  to 

equations  (3.2).  Equations  (3.1)  and  (3.3)  may  also  be 
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expressed,  in  terms  of  amplitude  and  phase  angle,  as 


u(x)  =  ajo)  + 

/  H^(n)  s/n  n(A  -  (x) 

n  - 1 

(3.4  a) 

v(a)  ^  a^(o)  + 

^  sin  n(A  -  0-) 

n  it 

(3.4  b) 

with  the  amplitude  and  phase  angle  of  equation  (3.4a)  given 
by 


b^{A) 


and  Veerin')  -  tan 


bjn) 


(3.5) 


Now  the  copower  series  of  U  and  V  may  be  defined, 
after  Kahn  [5] , 


^  I  U(^) 


a  (0)  Q  (0)  + 

U  V 


U-vr 


(n) 


(3.6) 


where  it  may  be  verified  upon  combining  equations  (3.4) 
and  (3.6),  that  the  copower  spectrum,  or  cospectrum,  has 
the  form 


Cjb)  =  iHjn)H>)  COS  n(cr(n)  -  o;(n))  (3.7) 

Here  the  symbols  H(n)  ,  and  (J(n)  ,  are  defined  as  in  equa¬ 
tions  (3.5).  It  may  also  be  noted  that  for  the  latitudinal- 
mean  values  of  U(A)  or  V"  (a)  ,  the  cospectrum  L-uv  in 
equation  (3.6),  becomes  simply  the  power  spectrum  of  u\a) 
or  v\a) 

z 

iH(h)  or  FJ.ih)  -  iHin)  (3.8) 

Now  the  value  of  C  uv  at  any  given  wave  number  may  be  shown 
to  be  independent  of  its  values  at  any  other  wave  number. 
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Equations  (3.7)  and  (3.8)  may  therefore  be  used  to  ob¬ 
tain  the  contribution  of  a  specific  vv/ave  number  to 
either  the  cospectrum  or  the  power  spectrum,  as  the  case 
may  be. 

A  useful  expression  for  may  be  obtained  by 

subtracting  ^4^^^  from  the  left  side  of  equation 

(3.6),  after  which 


U  V  = 


Jh) 


Making  this  substitution  into  the  transfer  equation  (2.2) 
gives 


T 


c) 

do 


U  \ 

ccr^j 


(3.3) 


Thus  T’  may  be  expressed  as  a  sum  of  contributions  for 
each  wave  number. 

A  more  convenient  form  of  equation  (3.9)  will  now  be 
derived.  Making  use  of  the  geostrophic  approximation  to 
the  horizontal  flow,  the  eastward  and  northward  compo¬ 
nents  are 


U 


3 


3 

fa  c) 


V  ^ 


9 

facoscj) 


(3.10) 


Here  e  is  the  acceleration  of  gravity,  f  the  Coriolis 
parameter,  and  2  contour  height.  The  other  symbols  are 
as  previously  defined.  Denoting  the  n-th  harmonic  of 
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the  height  field  at  latitude  (j)  as 

2{n)  =  A(n)  COS  n(A-H^(n'))  i3.lly 

where  A(h)  and  Vln)  are  the  amplitude  and  phase,  respec¬ 
tively,  at  latitude  (])  ,  equations  (3.10)  are  expressible 


as 


U  Co)  -  ~  Y~ 
fa. 


COS  n ( "A  -  Hin))  +  sin  n (^ - v3. 1 2) 

^  y  c)  d)  ^ 


V(n)  ^  sin  n(A-'f(n)) 

TQCOS  4) 


v3.i3y 


Henceforth  the  wave  number  h  will  be  understood  in  A'n)  and 
Equation  (3.12)  may  also  be  written  in  the  form 


u,w  = 


lA 

y 

where  the  angle  ~X_  is  given  by 


sin 


,,0  •  1*1  y 


X  -  t 


an 


/  A 


V 


ad/  ao 


i,3M5) 


and  is  therefore  the  phase  difference  between  and 

.  Using  the  expressions  (3.13),  (3.14),  (3.15)  and 

(3.7),  equation  (3.9)  takes  the  form 

COS  0  3  (  Cl 


Tin)  =-1, 


[XtlSl 

i 

\f\o6  0  / 

m 


a  c)d)Vcosdi 

Equation  (3.16)  will  be  termed  the  energy  transfer  spec¬ 
trum.  By  simplifying  the  quantity  within  the  square 
brackets,  the  contribution  for  a  single  wave  number 
becomes 

‘AX  ( I^  M 


i^3. 16) 


1  ^ 

(n) 


y  c 


n 

iaf  ‘cosdv 


III) 


5  h-) 
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4.  The  24-hour  change  of  eddy  kinetic  energy. 


An  approximation  to  the  change  of  eddy  kinetic  energy 
in  a  closed  system  resulting  from  all  of  the  processes  of 
equation  (2.1)  is  readily  made.  The  eddy  kinetic  energy 
per  unit  mass  at  any  position  is  given  by 


KE(n)=  Y(n)  =  ^a(p)  + 


(4.1) 


Using  equations  (3.8),  (3.12)  and  (3.13),  one  obtains  the 

latitudinally-averaged  or  spectral  values  at  latitude  cj) 


U 


^(n)  =  R(h)  = 


\-(n)  = 
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(4.2) 


Substitution  of  equations  (4.2)  into  (4.1)  yields  the 
zonally-averaged  geostrophic  eddy  kinetic  energy  per  unit 
mass.  For  a  closed  system  of  mass  /A, 


K  Em  =  Y 


+  i-^f  s 


(4.3) 


30 1  '  \n/T  30  j  \to5  0 

The  mean  24-hour  change  is  found  using  the  simple  differ¬ 
ence,  shov?n  symbolically. 


AKL(r\)  -  KE(h'),  —  Kl(im  .  (4.4) 

do.-/  Ui  '■ 

5.  Finite-difference  forms. 

The  element  of  mass  S m  of  equations  (3.17)  and  (4.3) 
may  be  expressed  initially  as 

5m  -  p  cx.^  cos  <t>  S(p  Sz  (5.1) 
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which,  upon  substitution  from  the  hydrostatic  equation, 
becomes 


<5,.  =  5x54- S-P  (5.2) 


where  o-p  ,  the  pressure  increment  is  considered  positive. 

In  a  series  of  papers,  Saltzman  [10],  and  Saltzman  and 
Fleisher  [11,  12]  have  been  careful  to  simulate  a  closed 
system  by  using  a  large  portion  of  the  northern  hemisphere 
at  500  mb,  for  example,  from  15°  N  to  80°  N.  In  this  thesis, 
it  seemed  advantageous,  from  the  point  of  view  of  a  possible 
forecast  tool,  to  investigate  kinetic-energy  changes  in  ten- 
degree  latitude  belts,  even  though  such  belts  may  not  con¬ 
stitute  closed  systems. 

The  major  purpose  of  this  study  is  to  investigate  pos¬ 
sible  relationships  between  T(T.)  and  ziKE(n)at  the  500-mb 
level.  Combining  equations  (3.17)  and  (5.2)  leads  to 

2fr 


3 


+  u 


tan  (t) 


t-adi  f,  o 

Integration  of  (5.3)  with  respect  to  longitude  will  increase 
the  constant  term  by  a  factor  of  277  .  The  remaining  inte¬ 
grations  are  with  respect  to  pressure,  where 


l-pz-  P..  -  p 

centered  at  500  mb,  and  with  respect  to  latitude,  where 
2  A(])  =  10  degrees,  or  TT/lS  radians. 

Employing  finite  differences  and  the  geostrophic 
approximation,  the  quantity  ^  in  equation  (5.3)  becomes 
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D4> 


3 


2  rLa^<t) 


'l  _ 

l^t>  As  \ 

VSin((P  1 1  ‘^l 

\  S)  h  '  A . 

A(I> 


(5.4) 


where  Z  is  the  height  of  the  mean  height  contour  at  the 
latitude  indicated  by  the  subscripts.  Moreover  U  ,  at  the 
central  latitude  (P  ,  may  be  expressed  as  the  arithmetic 
mean 


U  -  - 


/  ^ 

-h 

1  £<t  ™ 

Vs)  n  (d  -  4  <)/ 

2 

- 

(5.S) 


Substituting  (5.4)  and  (5.5)  into  the  integrated  form  of 
equation  (5.3),  the  finite  difference  form  may  be  stated  as 


T(n)  =  4t  <3^ n  A 


-a  1  -  ‘'(j) 


_ 

.  Sm  ((t)t2  5y  V  -  ‘  n  (^r  -  4 
’  4  4)  J 


+ 


<l).f  ^(t  \  -j-  “'*■ 


Sin(a.wV!.5)/  ^Sin(i 

L  ~  a 


-  4  .^y 


Combining  equations  (4.3)  and  (5.2)  yields 


KE(h)  - 


4  (^i^)  J 


nA|t)  ->  edJ 


Sa  Sfbct  {S.l) 


which,  when  integrated  as  before,  takes  the  finite-difference 


form 


KE(n)  - 


\  2ACl)  )  ' 


M 

vcosd). 


2  A  (t 

Za^ 


(.5.S) 


where 


A  '^A^  +  Adi-sj  J 


ahd  2  Au¬ 


lt  should  be  noted  that  the  same  integrations  have  been 
performed  in  arriving  at  equation  (5.8)  as  for  (5.6).  The 
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units  of  are  therefore  those  of  kinetic  energy  per 

mb  and  per  10°  latitude  belto 
6.  Data  processing  and  computation. 

Data  for  the  evaluation  of  equations  (5.6)  and  (5.8) 
were  taken  from  a  Fourier  analysis  of  the  500-mb  pressure 
surface  which  was  made  by  the  Naval  Weather  Research 
Facility  [13]  <,  The  Facility  used  the  Historical  Series 
500-mb  hemispheric  charts  for  the  winter  season  of  1947- 
48 o  Mean  500-mb  heights,  and  the  amplitude  and  phase  for 
each  of  the  first  18  wave  numbers  were  computed  Values 
were  available  for  five-degree  intervals  from  15  degrees 
to  70  degrees  North  latitude.  Data  selected  for  this 
study  were  for  the  33-day  period  from  31  December  through 
1  February. 

Because  of  time  limitations,  only  waves  1  through  12 
were  considered «  However,  results  of  Duggan  [3]  and 
Melhorn  and  Le  Dew  [7]  indicate  that  contributions  from 
wave  numbers  higher  than  ten  may  be  expected  to  be  rela¬ 
tively  small. 

In  the  case  of  phase-angle  data,  further  refinement 
was  necessaryc  Phase  angles  were  given  in  terms  of  dis¬ 
tance  from  the  prime  meridian  to  the  ridge  line  of  the 
wave  in  question.  Because  of  the  method  of  tabulation  of 
phase  angles  used  in  the  original  harmonic  analysis,  this 
value  did  not  always  pertain  to  the  ridge  nearest  the 
meridiano  Thus,  if  the  phase  difference  obtained  from  the 
reported  values  was  in  excess  of  one-half  wave  length,  it 
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was  necessary  to  subtract  one  wave  length  from  the  orig¬ 
inal  difference.  Fig.  1  schematically  depicts  this  process 
of  identification. 

Computation  of  equations  (5.6)  and  (5.8)  was  performed 
for  each  of  five  ten-degree  latitude  belts,  centered  re¬ 
spectively  at  65,  55,  45,  35,  and  25  degrees.  The  Control 
Data  Corporation  CDC  1604  computer  was  employed  in  these 
evaluations,  being  ideally  suited  because  of  its  high  speed 
and  large  capacity. 

7.  Results  and  conclusions. 

It  was  indicated  previously  that  if  integration  of  the 
transfer  equation  (2.2)  is  performed  over  the  entire  atmos¬ 
phere  ,  an  approximate  measure  of  the  transfer  of  kinetic 
energy  between  the  eddies  and  the  mean  zonal  wind  is  ob¬ 
tained.  From  the  evaluation  of  this  integral  over  a  10- 
degree  latitude  belt,  as  is  done  in  this  study,  only  a 
partial  verification  of  this  approximation  can  be  inferred. 
Fig.  2  shows  the  33-day  mean-energy  transfer  spectrum 
'"r  (n)  at  the  central  latitude  of  each  10  degree  latitude 
belt.  A  positive  value  of  Hn) indicates  a  contribution  to 
the  kinetic  energy  of  the  eddy  flow  from  the  mean  wind , 
whereas  a  negative  value  indicates  that  the  eddies  supply 
kinetic  energy  to  the  mean  wind. 

Perhaps  the  most  striking  feature  of  the  spectra 
shown  in  Fig.  2  is  the  change  in  sign  of  the  transfer 
between  35°N  and  45°N.  At  the  lower  latitudes,  transfer 
from  the  eddies  to  the  mean  flow  is  indicated.  This  is 
in  keeping  with  the  findings  of  Kao  [6]  who  shows  a  maximum 
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Figure  1.  Schematic  il'lustration  of  the  identification 
•  of^ corresponding  wave- pairs  at  neighboring 
latitudes.  •  . 
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convergence  of  northward  transport  of  kinetic  energy  at 
500  mb  up  to  approximately  36^N,  which  is  in  the  vicinity 
of  the  zone  of  maximum  westerlies.  Thus,  at  latitude 
25  N  and  35*^N  the  transfer  spectrum  appears  to  assist  in 
the  maintenance  of  mean  zonal  winds.  It  should  be  noted 
that  waves  2  and  3  are  the  largest  individual  contributors 
at  these  latitudes. 

At  45°N,  a  change  in  the  sense  of  energy  transfer  may 
be  seen.  All  waves,  with  the  exceptions  of  7  and  11,  now 
receive  energy  from  the  mean  flow.  The  transfer  from  mean 
to  eddy  flow  occurs  at  55°N  and  65*^N,  although  at  65*^N  the 
magnitude  of  the  exchange  is  greatly  reduced. 

The  spectral  estimate  for  the  full  latitude  range  20*^N 
to  70°N  is  shown  in  Fig.  3.  In  this  instance  the  sign  of 
the  energy  transfer  spectrum  has  been  changed  to  -Id 
to  facilitate  comparison  with  the  results  of  Saltzman  [10] . 

In  Fig.  3,  a  positive  sign  now  indicates  a  contribution  to 
the  mean  wind  by  the  eddies. 

Results  of  this  study  agree  roughly  with  those  of 
Saltzman  [10] .  The  major  difference  between  the  two 
spectra  lies  in  the  net  transfer  for  all  of  the  waves. 
Saltzman' s  results  show  an  overall  transfer  from  the  eddies 
to  the  mean  flow,  ^ i.e. , .  ^  ,  )  >  0  ^  ,  while  in  the 

present  study  there  is  a  net  transfer  of  kinetic  energy 

r*  ^  ^  * 

from  the  mean  to  the  eddy  flow,  |  /  '  '  i  ^  ; .  This 

is  not  unexpected  in  the  light  of  Saltzman  and  Fleisher's 
[11]  quoted  values  of  the  standard  deviation  of  the  transfer, 
which,  in  general,  exceed  the  value  of  the  spectral  estimate. 
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It  is  reasonable  to  suppose  that  changes  in  at 

a  given  wave  number  may  be  accompanied  by  changes  in  the 
kinetic  energy  associated  with  other  wave  numbers  at  the 
same  latitude  and  at  adjacent  latitudes.  Accordingly, 

' — I — I  ¥ 

values  of  i  in)  at  latitude  45  degrees  were  linearly  cor¬ 
related  with  those  of  ^  -  '  for  wave  numbers  n  =  1, 

r  yn  t 

2,... 12.  In  addition,  I  (n)  was  correlated  with  values  of 
at  adjacent  wave  numbers  nil  and  n  2.  Corre¬ 
sponding  values  of  at  latitudes  35  and  55  degrees 

were  also  correlated  with  the  value  of  I  ^  at  45  degrees. 
Table  1  gives  a  schematic  representation  of  this  correla¬ 
tion  system. 


Lat 

55° 

.  KE(n-2) 

aKE  (n-1) 

AKE(n) 

/..KE(n  1) 

KE(n 

2) 

Lat 

0 

in 

KE(n-2) 

KE (n-1) 

.  KE (n) 

KE (n  1 ) 

KE  (n 

2) 

Lat 

35° 

AKE(n-2) 

aKE  (n-1) 

aKE  (n) 

.KE(n  1) 

A  KE  (n 

2) 

Table  1.  15-way  correlation  scheme, 
T*(n)  is  correlated  with  each 
aK.E.  shown. 


Table  2  shows  the  results  of  the  15-way  correlation 

_  ^ 

1  TV 

system  in  which  I  (n  at  45  degrees  is  correlated  with  each 
of  the  15  sets  of  values  of  AKE (n , nil , ni2 )  at  latitude  35, 

45,  and  55  degrees.  The  values  of  AKE  are  centered  12 

►  -  ^ 

hours  prior  to  the  time  of  '  .  Tables  3,  4,  and  5  show 
the  results  with  t  centered  12,  36,  and  60  hours  respec- 
tively  after  the  time  of  T"  n  .  The  tables  are  arranged  in 
the  schematic  fashion  shown  in  Table  1. 
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An  excellent  statistical  routine,  written  for  the 
CDC  1604  computer  by  Campbell  [2] ,  was  employed  in  the 
computation  of  the  correlations. 

In  proceeding  to  the  correlation  analysis,  it  was 
necessary  to  establish  confidence  limits  for  the  correla¬ 
tion  coefficient.  The  number  of  independent  pairs  in  all 
the  correlations  was  30.  According  to  the  z-test,  as 
given  by  Panofsky  and  Brier  [8] ,  a  correlation  coefficient 
of  .30  is  significant  at  the  90%  level  of  belief. 

The  physical  principle  to  be  tested  is  the  conversion 
of  the  transfer-energy  by  wave  number  n  into  the  kinetic 
energy  of  wave  number  n,  nil,  nt2.  Saltzman  [9,  equa¬ 
tion  (49)]  has  derived  a  transform  of  equation  (2.1).  His 
result  includes  all  of  the  terms  which  appear  in  (2.1),  in 
spectral  form,  and  in  addition,  terms  representing  inter¬ 
action  of  kinetic  energy  from  wave  number  n  to  all  other 
waves.  In  a  later  paper,  Saltzman  and  Fleisher  [12]  have 
estimated  the  rate  of  energy  interaction  at  500  mb.  They 
have  shown  that  energy  transitions  tend  to  occur  between 
a  middle  group  of  waves,  (n  =  6  through  10)  called  "cy¬ 
clone  waves,"  into  the  longer  wave  group  (n  =  1  through  5) 
and  the  short  wave  group  (n  =  11  through  15) .  The  present 
study  is  more  limited  in  scope,  in  that  the  range  of  waves 
considered  is  only  (n  -  2)  through  (n  t  2).  On  the  other 
hand,  investigation  by  means  of  correlations  between  speci¬ 
fic  wave  numbers  may  make  available  more  detailed  criteria 
for  the  energy  exchange  between  waves  at  latitudes  35°, 

45°,  and  55°/  in  terms  of  at  45  degrees  latitude. 
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For  the  purpose  of  such  an  investigation  correlation 
coefficients  from  a  particular  row  and  column  of  each  of 
Tables  2,  3,  4,  and  5  were  selected  in  accordance  with  the 
criteria  that  each  set  of  four  successive  coefficients 
contain  at  least ; 

(a)  two  significant  correlations 
or 

(b)  one  significant  correlation  together  with 
a  24-hour  or  48-hour  difference  in  corre¬ 
lation  of  ±.60. 

All  together  there  were  47  sets  of  successive  correlations 
in  Tables  2  through  5  which  satisfied  criteria  (a)  or  (b) . 

As  an  example,  if  the  row  and  column  just  mentioned 
correspond  to  (11,  9),  each  of  the  four  tables  would  be 
entered  with  the  central  n-value,  n  =  11,  and  column 
n  -  2  =  9.  This  particular  correlation  coefficient  would 
be  symbolized  IT  c^  .  If  a  row  corresponding  to  the  higher 
latitude  (55°N)  is  selected,  a  bar  is  placed  over  the  T 
(e.g. ,  F  ),  whereas  if  the  lowest  latitude  is  used,  a  bar 
is  placed  below  the  T  (e.g.,  F  ).  For  a  choice  of  a 
central  latitude,  no  bar  is  employed. 

For  example,  signifies  the  correlation  between 

t — r  T  ^  Q 

I  (II)  computed  at  45  N  versus  aKE(9)  computed  at  latitude 
35  degrees.  Of  this  particular  set  of  four  correlations, 
three  were  significant  at  the  90%  level  or  higher.  The  set 
of  four  is  plotted  in  Fig.  4a.  F"!  |  u  shows  a  similar 

pattern,  but  in  this  case  has  a  significant  positive  cor¬ 
relation  at  time  t  -  12  hours.  This  initially  positive 
significant  correlation  tends  to  confirm  the  hypothesis 
that  an  above-normal  value  of  T(  n )  is  associated  with  an 
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Flgiirb  4,  i  SotJuGntiO-l  cortolo-tion  coefficient  setc  of  T''(n)  vjith  AKE 
.jiijji  for  selected  central  n-values.  Times  are  in  hours  from  the  com- 
tlui|!:l!|  putation  time  t  of  T*(n)*  YaluoB,oflH>  .JO  are  significant  at 
tho  pO/S  level  of  belief. 
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above-normal  value  in  AKE(n)  at  the  same  wave  number. 

At  time  t  4-  12  hours,  however,  the  value  is  negative,  os¬ 
cillating  to  positive  and  negative  again  on  the  succeed¬ 
ing  two  days.  These  results  are  consistent  for  all  cor¬ 
relation  sets  meeting  the  criteria  (a)  or  (b)  in  n  >  10 
and  suggest  that  above-normal  value  of  "Tcn")  tends  to  be 
followed  by  below-normal  change  in  perturbation  kinetic 
energy  twelve  hours  later.  Several  of  the  correlation 
sets  examined  in  this  group  displayed  this  two-day  peri¬ 
odicity.  Certain  of  the  other  sets  in  the  short-wave 
group  lose  significance  at  time  t  +  36  and  t  -t  60  hours 
even  though  they  may  have  met  the  criteria  (a)  or  (b) . 

At  the  opposite  end  of  the  wave  number  scale,  there 
is  again  the  tendency  for  aKE  (2  or  3)  to  be  above-normal 
with  a  high  value  of  I  i  ,  as  shown  in  Fig.  4b.  Again 
ziKE  (2  or  3)  becomes  negative,  but  in  this  instance,  two 
days  later,  indicating  the  periodicity  of  the  kinetic 
energy  changes  at  waves  2  and  3  is  apparently  of  the  order 
of  4  days. 

At  wave  number  5,  there  is  a  marked  tendency  for  a 
•\u  (5)  to  be  negative  at  t  ■+-  12  hours.  A  tendency  for 
z-  (6)  to  be  negative  at  this  same  time  also  exists  as 

shown  in  Fig.  4c.  This  suggests  a  positive  correlation 
between  the  two.  However,  the  exact  nature  of  any  such 
relationship  is  difficult  to  determine  by  the  technique 
employed  here. 

At  intermediate  wave  numbers,  as  represented  by  n  =  8 , 

_  tends  to  be  positive  at  t  -t  12  hours.  At  the  same 

n  O 
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are  of  opposite  sign  as 


time,  however,  Tl  ^  ,  and  Y~ 

8,6  ' e ,  1 0 

may  be  seen  in  Fig.  4d .  This  result  points  to  a  general 
agreement  with  the  results  of  Saltzman  and  Fleisher  [12] , 
namely^  when  a  decrease  of  kinetic  energy  occurs  in  the 
cyclone  waves,  kinetic  energy  tends  to  be  exchanged  both 
upward  and  downward  in  the  wave  number  scale.  A  schematic 
illustration  of  this  exchanged  energy,  after  Saltzman  and 
Fleisher  [12]^  is  shown  in  Fig.  5. 


I  -  5 


; 


h  -  6  -10 


11'  h  -  It) 


Figure  5.  Estimate  of  the  kinetic  energy 
exchange  between  cyclone  waves  and 
longer  and  shorter  waves.  After 
Saltzman  and  Fleisher  [12]  . 


A  direct  correlation  analysis  of  the  kinetic  energy 
changes  between  the  various  waves  would  have  yielded  seven 
times  as  many  correlations  as  the  present  scheme,  so  that 
analysis  of  the  results  would  be  extremely  involved.  Also 
the  results  would  have  had  no  relationship  to  T* (n) ,  which 
function  affords  the  possibility  of  being  a  predictor 
based  on  physical  considerations. 

It  is  necessary  to  comment  on  the  relatively  small 
number  of  significant  correlation  coefficients.  For 
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example,  in  Tables  2  and  3,  there  occurred  27  and  30  sig¬ 
nificant  values,  respectively,  where  16  (or  10%)  could 
have  arisen  due  to  chance.  Some  of  the  reasons  that  a 
more  conclusive  number  of  significant  correlations  did  not 
appear  may  be  the  following: 

(a)  Of  the  several  macroscopic  eddy  stress 
terms  in  equation  (2.1)  only  the  T*  term 
was  considered. 

(b)  Twenty-four  hour  changes  of  KE  may  not 
behave  linearly  when  correlated  with  T*. 

(c)  In  many  cases  the  initial  correlation  be¬ 
tween  T* (n)  and  AKE(n)  at  t  -  12  hours  was 
negative  but  became  significantly  positive  a 
day  later.  This  is  indicative  of  the  im¬ 
portance  of  the  process  of  exchange  of  kin¬ 
etic  energy  between  the  waves  in  the  manner 
described  by  Saltzman  [12] .  Such  exchanges 
undoubtedly  contribute  to  the  initial  state 
of  the  kinetic  energy  spectrum. 

(d)  There  is  some  "noise"  present  in  all 
Fourier  analyses  of  meteorological  parameters. 
There  may  also  have  been  some  undetectable 
reproduction  errors  in  the  copy  of  the 
original  data  used  in  this  study. 

In  spite  of  the  shortcomings  of  the  linear  correla¬ 
tion  technique  employed  in  this  study,  it  is  felt  that  iso¬ 
lating  "significant"  correlation  sequences  of  the  type 


23 


shown  in  Fig,  4,  led  to  several  valuable  results.  Among 
these  are  : 

(a)  the  short  period  of  kinetic  energy 
changes  at  the  high  wave  numbers; 

(b)  the  longer  period  changes  at  lower 
wave  numbers; 

(c)  the  tendency  of  wave  number  8  (in  the 
center  of  the  cyclone  wave  group)  to  ex¬ 
change  kinetic  energy  with  both  waves  6 
and  10. 
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APPENDIX 


The  following  tables  list  the  results  of  correlating 

the  energy  transfer  function  T* (n)  computed  at  45°N,  with 

the  24-hour  kinetic  energy  change  AKE  computed  at  35°, 
o  o 

45  /  and  55  N.  In  addition,  AKE  is  varied  with  respect  to 
wave  number  and  time  as  indicated  by  the  respective  table 
captions.  Significant  correlation  coefficients  are  de¬ 
noted  by  an  asterisk. 
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Central 

•"n" 

n  -  2 

n  -  1 

n. 

n  +  1 

n  +  2 

1 

r-i^wm 

.02 

-.01 

.10 

-.14 

•-.05 

-.05 

.50* 

.01 

.57* 

-.51* 

.26 

.06 

.10 

2 

-.01 

-.21 

.05 

-.19 

-.10 

.06 

.05 

.00 

.10 

.20' 

-.10 

.08 

-.01 

5 

-.01 

-.24 

.25 

.07 

.59* 

.07  . 

-.05  • 

.57* 

.08 

-.09 

-.29 

.04 

-.05 

-.25 

.10 

•  4 

-.59*' 

.02 

.05 

.27 

.05 

* 

-.05 

.10, 

-.07 

-.21 

.  .01 

.15 

•^06 

.08 

'  .10 

.15 

5 

..'10 

-.02 

.28 

.07 

.07 

.00- 

.12 

.►05 

.08 

-.11 

-.15 

.5^ 

.05 

-.54* 

-.22 

6 

-.05 

.15 

.20 

-.02 

-.05 

.12 

.06’ 

.  .05 

-.56* 

-.15 

.21 

.26 

-.47* 

-.55* 

.04 

7 

-.08 

.06 

-.42* 

-.28 

.57* 

.15 

.00 

.00 

-.10 

-.01 

.40* 

-.15 

-.27 

.25 

.28 

8  ’ 

.5^*' 

-.22 

-.21 

.10 

.05 

.  .11 

.10 

-.25 

-.06 

-.11 

* 

.26 

CO 

o 

• 

1 

-.02 

.00 

-.56* 

.9 

.15' 

-.54* 

-.56* 

-.09 

.02' 

.22 

-.18 

-.15 

.07 

-.11 

*02 

.50* 

-.24 

-.10 

-.19 

10 

-.10 

-.54* 

.55* 

♦  20 

.21 

.02 

<M 

• 

-.08 

.18 

.08 

.12 

-.15 

-.09 

.18 

—  Tin 

11 

.11 

.08 

.50* 

'.04 

.  .01 

-.19 

.61* 

.48 

-.09 

-.15 

.12 

12 

.51* 

.15 

.07 

...... 

-.52* 

.44* 

-.05 

.:;Tahle'  ;2.  .  Correlations  ’  of  ,'T*  (n)  at.,45*^N  at  time  .t,  with  iCsKE 
. .  ..<  7  ■  at  n,  n  i  1,’  and'- li  ±  2,  a't  time  -t  -  -12  hours.  -  Sig 
•  '  '  pificant  . values  ^are 'raarX ad- ..With  an  Asterisk. 
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"n" 

n  -  2' 

n  -  1 

n 

n  +  1 

n  +  2 

•• 

.10 

.12 

-.52 

1 

.11 

.55* 

-.19 

-.05 
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.00 

-.04 

-.05 

-.'08 

-.06 

2 

-.08 

.53* 

-.29 

.04 

.53* 

-.08 

-.55* 

-.22 

-.08 

-.01 

.04 

-.57* 

-.08 

5 

-.04 

-.07' 

-.59* 

-.09 

-.28 

.52* 

.05 

-.24 

-.15 

-.08 

,  . 

-.02 

.00 

-.09 

.09 

-.14 

4  ■ 

-.14 

-.20 

-.06 

-.58* 

.  -.15 

.12 

.20 

.06 

.05 

-.05 

-.01 

.17 

.12 

-.54* 

-.06 

5 

.29 

-.07 

-.71* 

-.16 

-.19 

-.10 

.07 

.04 

-.02- 

-.11 

.17 

-.21 

—.10 

—.06 

-.05 

6 

.06 

-.02 

-.52* 

-.51* 

^.06 

..05 

-.10 

.27: 

-.16 

-.52* 

V. 

-.2a 

-.04 

-.05 

^.21 

-.02 

7 

.12 

-.07 

.11 

.53* 

.05 

-.21 

.15 

—.02' 

-.52* 

-.17 

-.21 

.15 

.20 

-.19 

-.62* 

8 

-.45* 

—.07 

.49 

.10 

-.08 

.06 

.19 

-.08 

.21 

-.15 

.06- 

.04 

-.45* 

-.07 

.29 

9 

.01 

.11 

.59* 

-.48* 

.02 

.10' 

.54* 

-.12 

-.12 

.05 

r  r‘ 

.08 

-.25 

.15 

-.11 

.19 

10 

-.54* 

-.07 

-.29' 

-.58* 

-.15 

.27 

-.19 

.27 

-.45* 

-.42* 

• 

-.14 

.02 

fTN 

CM 

• 

-.09 

— 

11 

.16 

-.53* 

-.44* 

-.14 

-.57* 

.25 

-.53* 

-.27 

— 

-.08 

.24 

-.17 

12 

-.15 

-.10 

.04 

..... 

.16 

-.52* 

..’21 

— , 

-Table-  3.  ’  Correlations  ‘of  T*  (n)  atUs^N-and  time  t,  with 
^(KE)  at  n,i,.  JL‘,  and/nji- 2^  at  tlme,^.rt;i2  hoyr 
‘  •'  ,  .  Significant,  values  are  marked  with  an.  asterisk 

'29-  ' 


Central 

"n* 

n  -  2 

n  -  1 

n 

n  +  1 

n  +  2  , 

-.10 

-.06 

-.21 

1 

— 

.19 

-.05 

.16 

— 

-.18 

-.25 

-.55* 

1 

-.24 

.00 

-.15 

2 

.21 

-.09 

-.16 

-.14 

- 

-.15 

-.08 

.05 

.00 

• 

-.27  , 

.07 

.14 

-.01 

5  .. 

.11 

.09 

.10 

-.18 

-.05 

-.05 

CVJ 

• 

1 

.25 

..55* 

.12 

.55* 

.00 

CM 

CM 

.06 

-.09 

4  . 

.29 

.18 

.04 

-.01 

■  .09 

.01 

.;i4 

.59* 

.04 

-.29 

.05 

.02 

-.20 

.18 

.16 

5 

.17/ 

.22 

.21 

.09 

.15 

.07 

.05 

-.21 

.05 

-.19 

-.21 

-.17 

■  -.04 

-.05 

.07 

6 

.  -.02 

-.14 

-.24 

-.15 

-.06 

-.22 

.15 

.50* 

.15 

-.12 

.05 

-.10 

-.48* 

.  .01 

.14 

7 

-.05 

-.25 

-.08 

.01 

.05 

.11 

.14 

-.04 

-.20 

.50* 

r.07 

-.44* 

.22 

.15 

.44* 

8 

'  -.08 

.14 

-.06 

-.06 

.00 

-.07 

-.15 

.12 

.07 

.,20 

.17 

-.05 

.41* 

.18 

-.01 

9 

.16 

.25 

-.02 

.58* 

.00 

.07 

-.04 

.11 

.06 

.15 

• 

-.29 

-.05 

-.25 

.11 

.05 

10 

.08 

.19 

.06 

.22 

.05 

-.00 

.09 

-.52* 

.07 

.19 

.18 

.08 

-.08 

-.18 

11 

-.09 

.55* 

.50* 

.17 

.41* 

-.14 

.28 

.02 

— - 

.14 

-.04 

-.09 

— 

— 

12 

-.11 

.11 

-.15 

-.05 

.08 

-.15 

— ■  — — 

Correlations  of  T’^(n)  at  45°N  ard  time  t  with 
.jKE  at  n  n  t  1  and  n  2/ at  time  t  36  hours. 
Significant  values  are  marked  with  an  asterisk. 


Table  4: 


Central 

•n«  • 

CM 

1 

n  -  1 

n 

n  +  1 

n  4*  2 

• 

ro 

0 

.05 

.59* 

1 

-.04 

-.28 

.07 

CM 

0 

• 

1 

.05 

.08 

•MW 

.06 

-.05 

.08 

.15 

2 

-.29 

.5^ 

-.20 

.17 

\ 

.10 

-.05 

-.09 

-.24 

0 

-.42* 

.24  . 

-.05- 

.15 

.06 

5 

-.06 

.15' 

.17 

.42* 

.05 

-.05 

.46* 

-.47* 

-.15 

.00 

,02 

-.06 

-.05 

.09 

.14 

4  • 

-.50* 

-.17 

.51* 

.05 

•  .15 

.56* 

-.04 

-.04 

.05 

.26 

-.01 

-.07 

.09 

.10 

-.04 

5 

-.57* 

.25 

-.05 

-.05 

.45* 

'  .18 

-.15 

.02 

.15 

.14 

.07 

-.09, 

.15 

-.05 

-.05 

6 

.04 

.06 

<r.11 

.11 

.05 

1 

• 

O 

.04 

-.41* 

.25 

.05 

-.02 

-.14. 

.41 

-.05 

.08 

7 

.02 

.04 

-.55* 

-.10 

-.45* 

.26 

1 

• 

O 

.05 

.02 

-.20 

.04 

.55* 

-.17 

.11 

.11 

8 

.21 

-.20 

.00 

-.45* 

-.09 

-.09 

-.26 

-.04 

-.14 

-.08 

, 

-.19 

.09 

.01 

-.05 

-.25 

%  ■ 

■  .04 

-.15 

-.18 

.05 

•  .15 

f 

-.27 

-.29 

-.01 

.10 

-.06 

..  .  “ 

.16 

-.22 

-.07 

.00 

.01 

.10  ^ 

-.11 

.05 

-.10 

-.25 

-.11 

— .  06 

.05 

.26 

.00 

.08 

* 

-.12 

CM 

O 

• 

1 

-.10 

.14 

11 

-.51* 

-.50* 

-.28 

-.07 

-.51* 

.18 

-.17 

.09 

.  V 

-.04 

.04 

.14 

T 

’  12 

-.12 

-.08 

.00 

.17 

.01 

.19 

Table  5.  Correlations  of  T’^(n)  at  45°N  and  time  t,  with 

/iKE' at  n,  n  +  .  l,  and  n  ±  2, at  time  t  +  60  h.ours. 
Significant  values  are  marked  with  an  asterisk. 
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